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Abstract. The negatively charged nitrogen- vacancy centre is a unique defect in 
diamond that possesses properties highly suited to many applications, including 
quantum information processing, quantum metrology, and biolabelling. Although 
the unique properties of the centre have been extensively documented and utilised, a 
detailed understanding of the physics of the centre has not yet been achieved. Indeed 
there persists a number of points of contention regarding the electronic structure of the 
centre, such as the ordering of the dark intermediate singlet states. Without a detailed 
model of the centre's electronic structure, the understanding of the system's unique 
dynamical properties can not effectively progress. In this work, the molecular model of 
the defect centre is fully developed to provide a self consistent model of the complete 
electronic structure of the centre. The application of the model to describe the effects 
of electric, magnetic and strain interactions, as well as the variation of the centre's fine 
structure with temperature, provides an invaluable tool to those studying the centre 
and a means to design future empirical and ab initio studies of this important defect. 



PACS numbers: 31.15.xh; 71.15.-m; 76.30.Mi 



Submitted to: New J. Phys. 



The electronic solution of the NV centre 



2 



1. Introduction 

The negatively charged nitrogen- vacancy (NV - ) centre in diamond is a promising system 
for many quantum information processing pQ, quantum metrology, and biolabelling 
applications |2J. These applications include secure quantum key distribution [3], 
quantum computing [I], Q-switching in solid state photonic cavities [5], magnetometry 
[6] , electric field sensing [7] and decoherence based imaging [8] . The significant interest in 
the centre is primarily due to its well documented capabilities of single-photon generation 
[9] , long-lived coherence [10] , spin coupling [11] and optical spin polarization and readout 
[12] . The observed properties of the centre include a strong optical zero phonon line 
(ZPL) at 1.945 eV [13j, an infrared ZPL at 1.190 eV [13] , a paramagnetic ground state 
triplet [15] . and a strain [IB] , Zeeman [17] and Stark [18] affected excited state triplet. 
Recent experimental studies have also provided new information regarding the centre's 
excited state fine structure, its temperature dependence [19] and the presence of dynamic 
Jahn- Teller effects [20J. 

The NV~ centre is a point defect of C^ v symmetry in diamond consisting of a 
substitutional nitrogen atom adjacent to a carbon vacancy (refer to figure [I]). The 
observable properties of the centre are consistent with a six electron model [21], where 
the electrons are postulated to consist of the five unpaired electrons of the nearest 
neighbour nitrogen and carbon atoms to the vacancy and an additional electron trapped 
at the centre. The even number of electrons yields an integer spin system and the use 
of spin resonance techniques [22J has confirmed that the electronic states are highly 
localized to the vacancy and its nearest neighbours and that the ground state is an 
3 A 2 triplet [15]. The high degree of localization supports the application of a molecular 
model of the electronic system of the centre, in which the centre's electronic states are 
described by configurations of molecular orbitals (MOs) initially constructed from linear 
combinations of the dangling sp 3 orbitals of the nearest neighbour carbon and nitrogen 
atoms using group theoretical arguments. 

The previous applications of the molecular model [2ll |23j |24] have successfully 
described the zero field splitting of the ground triplet state due to spin-spin interaction, 
the 3 E excited triplet state and its fine structure induced by both spin-orbit and spin- 
spin interactions, and many aspects of the interactions of the triplet states with electric, 
magnetic and strain perturbations. However, being a semi-empirical model, unless the 
molecular model is fully developed in order to reduce the model's parameters to the 
minimal set, the model has limited ability to make definitive predictions on aspects of 
the electronic structure which can not be directly observed and the design of a systematic 
method to measure the large set of parameters becomes ambiguous. As a result of the 
previous partially developed molecular models, there has been a number of continual 
points of contention regarding the electronic structure of the centre. In particular, the 
contention surrounding the energetic ordering and positioning of the dark l E and 1 A\ 
(and possibly l E') singlet states that are thought to exist between the ground and 
excited triplet states. 
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Figure 1. Schematics of the nitrogen- vacancy centre and lattice depicting: the 
vacancy (transparent); the nearest-neighbour carbon atoms to the vacancy (black); 
the substitutional nitrogen atom (brown); and, the next-to-nearest carbon neighbours 
to the vacancy (white) . The adopted coordinate system and carbon labels are depicted 
in the right schematic. 

The locations of the intermediate singlet states are critical to developing an 
understanding of the process of optical spin polarization [12], which is the principle 
property of the NV~ centre that underpins the majority of its important applications. In 
the process of optical spin polarization, the population which has been optically excited 
from the ground to the excited triplet state is believed to decay non-radiatively from 
the excited triplet via the intermediate singlets to the ground triplet state in such a way 
that the m s = ±1 sub- levels of the excited triplet state are preferentially depopulated 
and the m s = sub-level of the ground triplet is preferentially populated. Consequently, 
after a short period of optical excitation the centre becomes spin-polarized into m s = 
population. This paper does not aim to describe the spin polarization mechanism, but 
instead provide the detailed model of the coupling of the intermediate singlet states 
and the triplet states due to spin-orbit interactions, which will form the foundations for 
future studies of the spin polarization mechanism. 

There are a number of other properties of the NV - centre that require a fully 
developed model in order to be satisfactorily explained. These include the Stark effect 
in the ground state triplet [3 EH], the small anisotropy of the effective electronic g-factor 
of the ground state triplet [26] . the strain splitting of the infrared transition between the 
intermediate singlet states [14J, and the presence of the Jahn- Teller effect in the l E and 
3 E [20]. Each of these properties requires the detailed treatment of electronic Coulomb 
repulsion, spin-orbit and spin-spin interactions which act to couple the electronic states 
of the centre and allow these properties to exist. The coupling of electronic states implies 
that the calculation of the effects of electric, magnetic and strain perturbations must be 
conducted in the complete basis of the centre's electronic states. 

In this article, the molecular model of the NV~ centre will be fully developed to 
provide an electronic solution that is experimentally testable and offers explanations 
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for many of the remaining questions regarding the centre. The model will be based 
upon previous applications of the molecular model and utilize invaluable ab initio and 
empirical results to draw conclusions and identify parameters which are known. The 
electronic Coulomb, spin-orbit and spin-spin interactions will be treated to determine 
the energies and couplings of the electronic states. Spin-orbit and spin-spin interactions 
will be treated using perturbation theory in order to produce simple energy and coupling 
coefficient expressions in terms of the minimal set of model parameters. Each of the 
parameters are provided as one- and two-electron matrix elements of the electronic 
interactions, allowing unambiguous evaluation by future ab initio studies. Electric, 
magnetic and strain interactions are also treated in order to allow future experiments 
to be designed to measure the remaining unknown parameters. The treatment of 
the interactions will also provide the foundations to develop an understanding of spin 
polarization, the Jahn- Teller effect and the temperature dependence of the centre's 
properties. 

2. The orbital structure 

Adopting an adiabatic approximation and considering the nuclei of the crystal to be 
fixed at their equilibrium coordinates R corresponding to the ground electronic state, 
the electronic Hamiltonian of the NV~ centre can be defined as 

H NV = ^2^ + V Ne (ri, R ) + V so (Xi, R ) + Vee(xi, *i) + V ss (xi, Xj) 

i i>j 

(1) 

where Xj = (r*j, s*j) denotes the collective spatial and spin coordinates of the i th electron 
of the centre, Tj is the kinetic energy of the i th electron, V^e is the effective Coulomb 
potential of the interaction of the nuclei and lattice electrons with the electrons of the 
centre, V so is the electronic spin-orbit potential, V ee is the Coulomb repulsion potential 
of the electrons of the centre, and V ss is the electronic spin-spin potential. Note that 
nuclear hyperfine interactions have been ignored. As for most molecular and solid state 
systems, the first step in solving Hnv is to obtain the solutions of the one-electron 
Coulomb problem, 

h = f + V Ne (r,R ) (2) 

which will be the MOs of the centre. Using the MOs, a basis of many-electron 
configuration states that are solutions of ^ hi can be defined and the remaining one- 
and two-electron components of Hnv can be treated in this basis. 

At this stage, the C^v symmetry of the ground nuclear equilibrium coordinates can 
be employed to construct the MOs of the defect. Using the basis {n, ci, c 2 , c 3 } (refer 
to figure Q] for labels) of tetrahedrally coordinated sp 3 atomic orbitals of the nearest 
neighbour carbon and nitrogen atoms to the vacancy, the MOs can be constructed as 
linear combinations of the atomic orbitals (LCAOs) with definite orbital symmetry. This 
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procedure has been conducted by a number of authors and an example of the resulting 
set of MOs [29J is 

ai(N) = n, ai(C) = 1 (a + c 2 + c 3 - 3S nc n), 

e * = vviW 2cl ^ 2 ~ C3) ' e " = 7r^s (C2 ~ C3) (3) 

where S nc = (n\ci) and 5 CC = (c\\c2) are orbital overlap integrals. 

Clearly the LCAO method is a highly approximate method of constructing the MOs 
as it uses a restricted basis set and does not consider the interactions between the MOs 
of the defect centre and the electron orbitals of the remainder of the crystal. Therefore, 
the key objective of the LCAO method is not to produce an accurate description of the 
MOs, but to produce the correct number of MOs of a particular symmetry type and 
to estimate their energy ordering. The results of ab initio studies j27[ [3Q1 EI] can be 
used to confirm the number of MOs of each symmetry type and their energy ordering. 
The majority of ab initio studies agree that there exists three MOs {a\, e x , e y } within 
the bandgap of diamond and that these resemble the highly localized MOs of ([3]), with 
additional contributions from atomic orbitals at the next-to-nearest neighbours and 
beyond. Furthermore, the studies show that the ai(N) and a%{C) MOs have mixed to 
form di and a[ such that a[ has been forced downwards in energy into the diamond 
valence band and a\ has significant contributions from both the nitrogen and carbon 
atomic orbitals. 

Using the six electron model of the NV - centre, the a[ MO will be completely 
filled by two electrons in the ground afafe 2 and first excited afa\e 3 MO configurations. 
There are several second excited MO configurations: afe 4 , a[a 2 e 3 and a^aie 4 . Due to 
the estimated proximity of the a\ MO to the valence band [291 [31], the first two second 
excited MO configurations afe 4 and a^afe 3 could be close in energy. The states of 
these second excited MO configurations could mix with the states of the ground and 
first excited states and affect their energies, but are not expected to play a significant 
role in the centre's properties themselves. Consequently, only the three MOs within the 
band gap are expected to contribute to the observable properties of the centre and only 
the ground and first excited MO configurations will be treated in detail in this work. 
Schematics of the three MOs in the region of the vacancy and their energy ordering is 
depicted in figure [2j 

The configuration states of the ground and first excited MO configurations of the 
centre are constructed by first defining linear combinations of products of four MOs 
that transform as a particular row of an irreducible representation of the C% v group. 
Using the irreducible representations and Clebsch-Gordon coefficients defined in [23J, 
examples of the symmeterized linear combinations of products of two MOs are 

<j)A 1 { a i a i) = a i fl i' 4>E,x( a i- e ) = a i e z' 4>E, y (aie) = atfy, 
0Ai(ee) = -^={e x e x + e y e y ), 0A 2 (ee) = -^=(e x e y - e y e x ), 
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Figure 2. Schematics of the three NV centre molecular orbitals (MOs) responsible 
for the centre's observable properties in the vicinity of the vacancy and their energy 
ordering. Red and blue components represent positive and negative contributions to 
the MO respectively. 



1 -1 
(j>E,x( ee ) = -^( e x£x ~ e y e y ), (j>E, y (ee) = -j=(e x e y + e y e x ) (4) 

The symmeterized combinations of products of four MOs can be constructed by- 
repeating the process used to construct the above products of two MOs. Once the 
symmeterized products of occupied MOs corresponding to each configuration state 
are constructed, the configuration states are formed by performing a direct product 
with an associated spin state and transforming the result into a linear combination of 
Slater determinants. The resulting configuration states $ C j^-s,m a nave definite orbital 
symmetry (j denoting irreducible representation and k denoting row of the irreducible 
representation), total spin S, and spin projection m s , and are explicitly contained in 
table [TJ Note that the construction of the configuration states in this way is completely 
analogous to LS coupling in atomic structure, where the atomic states are constructed to 
have definite orbital (L, mi) and spin (S, m s ) quantum numbers prior to the introduction 
of spin-orbit interaction. The configuration states may also be expressed in terms of 
holes rather than electrons and these expressions are contained in table IA1I However, 
the hole formulism will not be used in the remainder of this article. 

The configuration states are solutions of ^ hi, with each of the states of a MO 
configuration having the same eigenenergy as depicted on the left hand side of figure |3j 
Employing the Wigner-Eckart theorem 



(<f>f,9\dp,q\<t>j,k) 
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V 






Q 





(h\\o P m (5) 



where O is a tensor operator, (g, q, k) denote the rows of the irreducible representations 
{f,p,j) of the Cs v group respectively and (||||) is the reduced density matrix element, 
the eigenenergies of each MO configuration can be expressed in terms of reduced density 
matrix elements involving the MOs 

ale 2 : 2(a 1 \\h\\a 1 ) + 2(e\\h\\e) , a x e 3 : (ai||^||ai) + 3(e||^||e) (6) 

The introduction of the Coulomb repulsion potential ^ i>:) Ke( x i, x j) splits the MO 
configurations into distinct triplet and singlet energy levels. The diagonal matrix 



The electronic solution of the NV centre 



7 



Table 1. The configuration and spin-orbit states of the NV - centre expressed in 
terms of Slater determinants of the molecular orbitals. Second quantization notation 
has been adopted to denote the occupation of the molecular orbitals in each Slater 
determinant in the order \aiaie x e x e y ey}, where an overbar denotes spin-down. 
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elements of the ground MO configuration triplet and singlets are 

3 A 2 : Co + (<f) A2 (ee)\\V ee \\<f) A2 (ee)) 
l E: Co + iMeeWeeWMee)} 
1 A 1 : Co + iMeeWeeWMee)} 

where 

C = (<pA 1 (a 1 a 1 )\\V ee \\(l) Al (ee)) + 4 : ((j} E {axe)\\V ee \\(j) E {a 1 e)) - 2{<f) E (a 1 e)\\V ee \\(f) E (ea 1 )) 

Hund's rules indicate that the 3 A 2 triplet is the ground electronic state, which 
implies that ((j) A2 {ee)\\V ee \\(j) A2 {ee)) < {(j) E {ee)\\V ee \\4) E {ee)), ((p Al (ee)\\V ee \\(f> Al (ee)}. The 
difference in the singlet diagonal matrix elements is 

e = {(j> Al {ee)\\V ee \\4> Al {ee)) - (^(ee)||K e ||^(ee)) 

= 2(e x e x \V ee \e y e y ) = 2 / p xy (f 1 )V ee (f 1 ,f 2 )p 3: y(f 2 )dr 1 df 2 (7) 



where Pxy{r) = e a;(r) e 2/(r)- The above integral is a standard exchange integral of a 
charge distribution p(r), which has been proven to be positive definite. One such proof 
[33] from electrostatics uses Green's identity to show that 

/ p(n) _ 1 ^r p{f 2 )df x df 2 = — I \E{r)\ 2 df> (8) 
J Fi — r 2 | 4tt J 
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where E(f) is the electric field generated by the finite charge distribution p(r). Thus, 
considering just the diagonal Coulomb matrix elements, the l Ai singlet must be higher 
in energy than the l E singlet. It has been shown that the difference in the diagonal 
Coulomb matrix elements of the X E singlet and 3 A 2 triplet (((f> E (ee)\\V ee \\(f> E (ee)) — 
(0a 2 ( ee ) 1 1 ^ee 1 1 0^4 2 ( ee ) ) ) 15 equal to the difference between the two singlets e [3lj, thereby 
confirming the ordering indicated by Hund's rules and implying that the states of the 
ground configuration are equally spaced prior to coupling with states of the higher MO 
configurations. In the first excited MO configuration, Hund's rules also indicate that the 
3 E triplet has a smaller repulsion energy than the l E' singlet and a similar argument as 
used for the singlet ordering in the ground MO configuration has been shown to confirm 
this ordering [35] . 

Since the only configuration states of the ground and first excited MO configurations 
that have both the same orbital symmetry and spin state are the l E and 1 E' singlets, 
these are the only states that are coupled by the Coulomb repulsion potential, and 
all of the other states are solutions of the orbital components of H^v, H Q = hi + 
Yli>j Ke( x i, x i)- The coupled E singlet states Q^k-fi.o can ^ e expressed as 

$ ^, fc; o,o = N K [$ c E , ik . 0>0 + K$ c Etmo ] (9) 

where k = x,y, the coupling coefficient k is a function of the Coulomb repulsion matrix 
element 

(W^IW = <<Maie)HK e ||<Mee)> (10) 

and N K = (1 + |k| 2 ) is the normalization constant. The interaction of the two E 
singlet states will also force the singlets apart in energy, shifting the lower 1 E singlet 
lower in energy towards the ground triplet state, and shifting the higher 1 E' singlet 
further higher in energy than the excited triplet state. 

Defining Ej-s to be the orbital energies of the configuration states, the derived 
orbital structure prior to the introduction of spin-orbit and spin-spin interactions is 
depicted in figure [31 The known optical ZPL (1.945 eV) and infrared ZPL (1.190 eV) 
transition energies are also included in the figure. Notably, the energy separations 
of the triplet and singlet states have not yet been observed. As indicated by (EJ), 
all of the configuration states of the ground MO configuration have the same nuclear 
equilibrium coordinates (ignoring the Coulomb coupling of l E and l E'). Likewise, all of 
the configuration states of the first excited MO configuration also have the same nuclear 
equilibrium coordinates, but these differ from those of the ground MO configuration. 
Defining the energy of the ground 3 A 2 state to be zero (-EU 2 ;i = 0)> the energy of the 
1 A 1 singlet can be expressed as E^o ~ E E . Q + 1.190 eV, where the infrared ZPL has 
been directly used since, correct to first-order in k, both 1 A\ and l E have the same 
nuclear equilibrium coordinates. The energy of the excited 3 E triplet in the nuclear 
equilibrium coordinates of the ground MO configuration is E E -\ = E$ + 1.945 ~ 2.180 
eV, where E$ ~ 0.235 eV is the Stokes shift of the optical transition |28j. Thus, the 
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configuration energies E E . and Ee'-o of the E singlet states, alongside their Coulomb 
coupling coefficient k, are the first unknown parameters of the molecular model. 



"'E 



1.190 eV 



S s +1.945 eV 



Figure 3. Energy level diagram of the orbital structure of the NV - centre. The 
molecular orbital (MO) configuration energies are depicted on the left hand side and 
the splittings into singlet and triplets due to the introduction of electronic Coulomb 
repulsion is depicted on the right hand side. The dashed lines indicate the locations 
of the E singlet states prior to their Coulomb coupling, e = 2(e x e x \V ee \e y e y ) is the 
difference in the Coulomb energies between the ground MO configuration states prior 
to the coupling of the E singlets. E$ ~ 0.235 eV is the Stokes shift of the optical ZPL. 



A recent ab initio study [31] concluded that the upper l E' singlet energetically 
crossed the 3 E triplet as the nearest neighbour carbon and nitrogen nuclei were 
symmetrically displaced. The proceeding argument clearly shows that this can only 
be the case if there existed a strong Coulomb repulsion interaction between the l E' 
singlet and a higher energy E symmetric singlet. Such an interaction would force the 
X E' singlet lower in energy and if it were large enough, it could potentially overcome the 
repulsion from the lower 1 E singlet and the difference in the Coulomb energies between 
the 1 E' singlet and the 3 E triplet. The higher energy E singlet that produces this effect 
is most likely the 1 E" singlet that arises from the second excited MO configuration 
a\a\e 3 . It is also possible that the higher energy E singlet arises from conduction band 
states, however, it would not be expected that a singlet formed from such delocalized 
states would interact so strongly with the highly localized states of the centre. The other 
relevant second excited MO configuration afe 4 forms an 1 A' 1 singlet. Consequently, the 
singlets of these two second excited MO configurations could potentially force the E 
and A\ singlets of the ground and first excited MO configurations down in energy by 
different amounts, resulting in the 1 Ai crossing one of the E singlets. The most recent ab 
initio studies [31J1 EI] indicate that the ordering of the ground MO configuration singlets 
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depicted in figure E] is correct and without further evidence that the second excited MO 
configurations states significantly shift the positions of the lower singlets, the simple 
orbital structure of figure [3] will be assumed. Due to the parametric formulation of 
the molecular model, the failure of this assumption will not effect the validity of the 
expressions derived in this analysis. 

Another recent ab initio study [30J indicates that the 1 A 1 singlet is close to and 
possibly higher in energy than the 3 E triplet. As the quantitative values of the Coulomb 
repulsion matrix elements are not known, the proceeding analysis can not offer a 
definitive conclusion on the ordering of the 1 Ai singlet in relation to the 3 E triplet. 
The current understanding of the mechanisms of spin polarization and optical dynamics 
[2U [36] appears to strongly indicate that a singlet state is close, but lower in energy than 
the 3 E triplet. Therefore, further ab initio work is required to definitively determine 
the ordering of the l Ai, 3 E and l E' states. 

Thus far, it has been shown how the orbital structure of the NV~ centre can be 
derived from the definition of the centre's MOs. The configuration states that have been 
obtained are solutions of the orbital components H of the centre's Hamiltonian and the 
orbital structure was determined up to the two unknown energies of the E symmetric 
singlets. The analysis of the Coulomb repulsion matrix elements has offered an insight 
into the contention surrounding the ordering of the intermediate singlet states. In the 
next section, the spin-orbit V so and electronic spin-spin V ss potentials will be treated as 
first-order perturbations to H Q and the fine structure and mixing of electronic spin of 
the centre will be determined. 



3. Fine structure and the mixing of electronic spin 

The construction of the configuration states in the previous section to have well defined 
orbital symmetry aided in simplifying the treatment of the orbital components H Q 
of the centre's Hamiltonian. Since spin-orbit and spin-spin interactions act on both 
the electronic orbital and spin coordinates, their treatment can be likewise greatly 
simplified by constructing linear combinations of the configuration states that have well 
defined spin-orbit symmetry. This can be conducted by first defining the combinations 
of the S=1,0 spin states (\S,m s )) that transform as particular rows of irreducible 
representations of the C 3v group [23] 

^ = 10,0), S A2 = \1,0), 

S E>X = -|(|1, 1) + |1, -1)), S E ,y = ^(|1, 1) - |1, -1)) (11) 

In an analogous method to that used in constructing the combinations of products of 
two MOs that had well defined symmetry, the above symmeterized spin states can be 
used in conjunction with the configuration states to construct the symmeterized spin- 
orbit states k contained in table [lj where n denotes the energy level, and j and k 
denote the irreducible representation and row that the spin-orbit state transforms as in 
spin-orbit space. 
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The spin-orbit and spin-spin interaction potentials are 



2m 2 c 2 

_ ^09ef-B \^ Si " Sj ■ rij){fjj ■ Sj) \- r . p. -> / 19 n 

~ 4vr^ 2^ |4f " 2. a * ' ^ • ^ l^J 

where pi and s*j = {h/2)&i are the momentum and spin operators of the i electron, 



iy ty nr* 

% 3 3 1 1 



x^x + y^y + ZijZ (x, y, z unit coordinate vectors), g e is the free electron 
g-factor, and A and D are rank one and two orbital tensor operators respectively. The 
components of the orbital tensor operators are 

A = -X E , y x + \ E ,xV + X A2 z 

— \Da-l — De,x,1 De, v ,1 —De,x,2 

D= De, v> i —\£ > a 1 +De,x,i —De, v ,2 

\ —De,x,2 —De, v ,2 Da 1 



(13) 



where 





3z 2 


167T \fij\ 3 






3(4 - y 


327r|fy 3 


\ r ij 5 




3ZijXij 


167r|r*jj 3 





2 ) ,, „2,,2 



u E , x ,i - l3 l5 , - 167r| ^ |3 

Using the spin-orbit states of table [1] and the above definitions of the tensor operators, 
the application of the Wigner-Eckart theorem allows the computation of the matrix 
representations of the spin-orbit and spin-spin potentials in terms of one- and two- 
electron reduced density matrix elements. The matrix representations are contained 
in tables |2] and [3] and the reduced density matrix element expressions are contained in 
table HI Note that as only the spin-orbit states associated with the triplet configuration 
states have non-zero spin-spin matrix elements, the matrix representation of the spin- 
spin potential is presented in the reduced basis of triplet spin-orbit states. In order 
to maintain clarity, only the upper triangle of the matrix representations have been 
presented. The lower half can be inferred using the hermitian property of the potentials 
and the fact that each of the reduced density matrix element expressions are real. 

In calculating the non-axial orbital components of the spin-orbit potential, it was 
concluded that the reduced density matrix element (e||As||e) must be zero and so it was 
not included in the matrix representation of table [2j This follows from the fact that as 
Xe, x is a purely imaginary hermitian orbital operator and its diagonal matrix elements 
in the MO basis are (e x \\E,x\ex) = ^75 ( e I I-VeI l e ) — ~( e y I A^le^), these diagonal matrix 
elements must vanish, implying that (e||A,g||e) = 0. In a previous application of the 
molecular model to the fine structure of the excited triplet state [H] this conclusion 
was not made and the coupling of <&f E k and $g° E k was incorrectly assigned as arising 
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Table 2. The matrix representation of the spin-orbit interaction potential in the basis 

/<B so <b so <f> so <B so <f> so <f> so <f> so (f> so (f> so <b so (b so i> so <b so <b so \ 

The spin-orbit parameters Ay and A_l are defined in table 3J The lower half of the 
matrix can be obtained using the hermitian property of the spin-orbit potential. 
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Table 3. The matrix representation of spin-spin interaction potential in the triplet 

basis {'i-n ••Kx ,• 'K-,: ,- n%, y , n\x, n° E ,y, nx^nuJ- The spm-spm 

parameters are defined in table |4l The lower half of the matrix can be obtained using 



the hermitian property of the spin-spin potential. 

-2D 1>Al V2D 1>E>2 

D Ml V2D 1>E>1 D ltE>2 

D Ml V2D 1>Etl D 1>E>2 

D 2Al V2D 2>E>2 

D 2Al V2D 2jEj2 

-2D 2Al 

-2D 2Al 



D 2jAl - 2D 2tEtl 

D 2Al + 2D 2 ^ E>1 



from the non-axial spin-orbit interaction matrix element — iy2{e\ \X E \ |e) instead of the 
correct spin-spin interaction matrix element V2D 2:Ej2 as contained in table |3j 

As the spin-orbit and spin-spin matrix elements are observed to be of the order 
of MHz-GHz [T8l IT6] and the energy separations of the singlet and triplet states are 
expected to be of the order of meV - eV (~ 10 2 — 10 5 GHz), it is appropriate to 
treat the spin-orbit and spin-spin potentials together as first order perturbations to H a . 
The configuration energies Ej. s and spin-orbit states are then the zero order 

energies and electronic states of the perturbation expansion. The application of first 
order perturbation theory yields the fine structure energies E n contained in table |5] and 
depicted in figure H] as well as the coupling coefficients contained in table IA21 The first 
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Table 4. The spin-orbit and spin-spin parameters of the molecular model expressed 
in terms of one- and two-electron reduced density matrix elements containing the 
molecular orbitals. Known values |16j and estimated order of magnitude of unknown 



values are 


contained in the right most column. 




Parameter 


Expression 


Value 


A,, 


-i(e||A^ 2 ||e) 


5.3 GHz 


x± 


=j=l\ax\\\ E \\e) 


- GHz 




2^ A2 (ee)||I? Al ||0 A2 (ee)) 


2.87/3 GHz 


D\_E,1 


A(4> E {ae)\\D EA \\4> A2 {ee)) 


- MHz 


Dl,E,2 


-4{Mae)\\D Et2 \\4> A2 (ee)) 


- MHz 




((f> E (ae)\\D Al \\(f> E (ae)) - ((f) E (ae)\\D Al \\(f> E (ea)) 


1.42/3 GHz 


D2,E,1 


-2((q> E (ae)\\D EA \\Mae)) ~ (Mae)\\D EA \\cj) E (ea))) 


1.55/2 GHz 


D2.E.2 


2{{ ( j )E {ae)\\D Ea U E {ae)) - {c)) E {ae)\\D Et2 \\(j> E {ea))) 


200/ 'y/2 MHz 



order corrected electronic states &n'j,k are defined in terms of the coupling coefficients 
by 



where N n is the normalization constant. The coefficient n contained in table K2l is the 
coefficient that arises from the coupling of the degenerate &l° E ^ k and §l° E k states of the 
3 E triplet due to spin-spin interaction and is expressed in terms of the spin-orbit and 
spin-spin parameters as 

2V2D 2 ,e,2 

V = -r (15) 

A,, - 3D 2Al + [(A,, - 3D 2Al ) 2 + SIZ^I 2 ] 3 

and Nf) = (1 + \n\ 2 ) 1//2 is the associated normalization constant. Note that the 

(2) 

coefficients s n ,m do not denote coefficients that are truly second order in the spin- 
orbit and spin-spin parameters, but rather contain the first order products of the 
Coulomb k and degenerate spin-spin n coefficients with the other spin-orbit and spin- 
spin coefficients. As k and n are expected to be orders of magnitude larger than the 

(2) 

other spin-orbit and spin-spin coefficients, the s n ,m coefficients may potentially be only 
slightly smaller than the Sn} m coefficients. 

The spin-orbit and spin-spin parameters that are contained in the first order 
energies of table [5]have been determined in previous studies [TBI |T6] by directly observing 
the fine structure of the ground and excited triplet states and modeling the variation 
of the excited triplet fine structure with the application of strain and electric fields 
at low temperatures. The values of these known parameters are contained alongside 
the unknown parameters in table HI The unknown parameters are all involved only 
in the first order couplings of the spin-orbit states of the singlets and triplets and do 
not contribute to the first order energies. Consequently, the effect of these unknown 
parameters can only be indirectly detected through the observation of interactions of 
the centre with electric, magnetic and strain perturbations that can not be explained 
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Table 5. The electronic energies correct to first order in spin-orbit and spin-spin 
interactions. The energies calculated using the known parameters of table U are 



provided 


in the 


right column. 
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Figure 4. Energy level diagrams of the fine structure of the 3 A 2 and 3 E triplets of 
the NV~ centre: (a) 3 A 2 fine structure independent of strain; (b) 3 E fine structure in 
the absence of strain; and, (c) 3 E fine structure in the limit of large non-axial strain 
Expressions of the energies E„ are contained in tableland A = A|| + E> 2 ^a x and 
Q = |[(Aii — SD 2t A 1 ) 2 + 8|D2,_e,2| 2 ] 1 ^ 2 - The strain interaction parameter c^bI^jJ is 
defined in 03. 



by the zero order spin-orbit states and the coupling coefficients arising from the known 
parameters. Examples of such interactions are the presence of the Stark effect in the 
ground state triplet [3 [25] and the anisotropy of the ground state effective g-factor 
[26] . Frustratingly, it is precisely these unknown parameters which are likely to govern 
the process of spin-polarization as they allow forbidden transitions between the triplet 
and singlet states, however, this will not be discussed further in this article. Using the 
known values of similar parameters, the estimated order of magnitude of the unknown 
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parameters are also contained in table HJ 



4. Electric, magnetic and strain interactions 

The Stark shift V$, Zeeman effect Vz and strain potentials that describe the centre's 
interaction with electric E, magnetic B and effective strain £ fields are 

V S = J2^-E, V z = ^fJ2^ + ^)-B, % = (16) 

i i i 

where d\ = efi is the electric dipole operator and U = x is the orbital magnetic 
moment operator. Note that the above effective expression for the strain potential (in 
which strain is treated as an effective local electric field) is derived using the group 
operator replacement theorem [32] and can be unambiguously interpreted for uniaxial 
strain, but some care is required in its interpretation for non-uniaxial strain. 

Each of these potentials contain operators that act on just the electronic orbital 
coordinates and can be written in terms of sums of orbital tensor operators with definite 
symmetry properties 

di = d E>x x + d E , y y + d Al z, U = -l E , y x + l E , x y + h 2 z (17) 

As the spin-orbit states constructed in the last section are simple linear combinations of 
configuration states that have both well defined orbital symmetry and spin, it is straight 
forward to apply the Wigner-Eckart theorem to calculate the matrix representations of 
each of the tensor orbital operators in terms of a set of reduced density matrix elements. 
Since the expressions of the reduced density matrix elements are identical for all orbital 
tensor operators of the same symmetry, the matrix representations of general orbital 
tensor operators O = Oi °f each symmetry are contained in tables IA3I and IA4I The 
matrix representation of the electron spin operator S = Yli ^ i s contained in table 
IA5I These matrix representations can be immediately applied to calculate the effects 
of the different electric, magnetic and strain interactions in the basis of the zero order 
spin-orbit states. 

For example, the matrix representation of the strain potential in the basis of the 
zero order spin-orbit states of the 3 E triplet {® s 5 ° EtX , ®i° E>y , $6°E,a;> ®6°E, y i *7°^ 2 ' $ mJ is 

= 

/ <4,Ai£* -d btE £ y -d b:E £ x \ 

db,AiCz db,E%x —db,Ei y 

d b:Al £ z + d biE ^ x -d b:E £ y 

-d btE £ y d b>Al £, z - d biE £ x 

—d b , E £, y d btE £ x d b ^ Al iz 

\ ~d bjE £ x -d bjE £ y d b>Al £ z ) 

(18) 

where d b>Al = (ai\\d Al \\ai) + 3(e\\d Al ||e) and d b>E = (e||d^||e). Previous studies [TBI [T6] 
have shown that the above matrix representation yields a strain variation of the 3 E 
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triplet fine structure that agrees excellently with observation. A theoretical plot of the 
effect of non-axial strain is depicted in figure |5] and it clearly shows that non-axial strain 
splits the 3 E fine structure into upper and lower branches. Due to the D 2 ,e,2 spin-spin 
parameter that couples the m s = and m s = ±1 configuration states, there exists two 
level anti-crossings in the lower branch as the m s = changes from being the highest 
energy state of the lower branch to the lowest energy state [18]. In the vicinity of these 
anti-crossings the states of the lower branch are all significant mixtures of the different 
spin sub- levels [18J. As there is no anti-crossing in the upper branch, the mixing of the 
spin sub-levels in the states of the upper branch does not significantly vary with strain 

HB]. 




5 10 15 20 



dbMA (GHz) 

Figure 5. Plot of the variation of the 3 E triplet fine structure with the magnitude 
of the applied non-axial strain The plot was produced using the known spin- 

orbit and spin-spin parameters of table |4] and the matrix representation of the strain 
interaction (jT5)l . 



In the limit of large strain the upper and lower branches become identical, with the 
same mixing and splitting of the spin sub-levels. This can be demonstrated by calculat- 
ing the matrix representations of the spin-orbit, spin-spin and strain potentials in the 
basis of the 3 E configuration states „, $| , X;U , $| ..i _i, $L ; i 0l ®e,v;1,i, ®e, v x-i} 
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V 



-db,E^,x 











—db,E^,x 



------ 

db,Eiy 








~db,E_ix 



db,Eiy 



db^y \ 

d b . E iy 

_ _0_ 0_ _ _ dbMy 

d b ,E^x 
d b , E £x 
d b>E £ x ) 



(19) 



Clearly, if a £ x strain was applied such that db^x was much larger than the spin-orbit 
and spin-spin parameters, then the influence of the matrix elements of the upper right 
and lower left blocks would become insignificant and the configuration states would 
separate into identical E x and E y orbital branches. The off diagonal spin-spin matrix 
elements in the diagonal blocks are responsible for the mixing and splitting of the spin 
sub-levels within each branch. The fine structure of the 3 E triplet in the high strain 
limit is also depicted in figure HJ 

Strain also effects the l E' singlet and, only after the Coulomb coupling of the E 
sing lets, the l E singlet as well at first order. The matrix representations of the strain 
potential correct to first order in k in each of the corresponding basis sets {$3 
and {$l EiX , $( E J are 



<T> C \ 

E,xi ^3,E,yf 




d a ,A^z + ^d a ^ix 
— 2nd a ^y 



— 2Kd a} E^,y 
d a ,A^z — 2,Kd a ,Eix 



db^iz — {db,E + 2Kda tE )ix 

(d b , E + 2nd ajE )£ y 



(d btE + 2Kda )E )iy 



(20) 

where d a ,Ai = 2(ai||d J 4 1 ||ai) + 2(e| I |e) and d a> E = \d E \ |e). The strain splitting of 
the 1 E' singlet has not yet been directly observed, however the strain splitting of the 
1 E singlet has been observed [TJ] in the strain dependence of the infrared transition at 
1.190 eV, thereby confirming the coupling of the E symmetric singlets. Investigations 
of this strain splitting and the information that it can provide about the mixing of the 
E singlets are on going. 

The operator matrix representations constructed in this section may also be used 
in conjunction with the first order corrected spin-orbit states presented in the previous 
section to investigate the appearance of interactions at first and higher orders in 
the spin-orbit and spin-spin parameters. As discussed, these interactions provide 
important information on the unknown parameters and should be the subject of future 
investigations. 



5. The room temperature electronic structure 

The dynamic Jahn- Teller effect has been observed in the 3 E triplet [20] . This effect arises 
from the vibronic coupling of the spin-orbit states of the triplet |38j. The presence 
of the dynamic Jahn- Teller has the consequence that with increasing temperature E 
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symmetric phonons drive transitions between the states of the triplet. The phonon rates 
are dependent on temperature and firstly are observed to give an unusual temperature 
dependence of the optical linewidths and at higher temperatures the phonon transition 
rates become much greater than the observable lifetimes of the states |39j. At this point, 
any population which is excited into one of the states is distributed to the other states 
that it is coupled with before any radiative decay can occur. If the phonon transitions 
distribute equal population to each of the coupled states then the radiative transition 
will have the average energy of the coupled states. This process has been previously 
discussed in the literature as orbital averaging [T9] . 

Since electron-phonon coupling is an orbital operator, the matrix representations 
ffT8]) and f fl9|) of an orbital operator in the 3 E triplet can be used to determine the 
selection rules of the phonon transitions. The allowed transitions are depicted in figure [6] 
for both the low and high strain cases. It is clear that for both the low and high strain 
cases, the phonon transitions will distribute population between the states of a particular 
spin projection, and since the transition matrix elements are identical for each transition, 
there will be approximately equal population distributed to each of the states of a given 
spin sub-level. Due to the coupling of the states with m s = ±1 spin projections by 
the D^e,i spin-spin parameter, the population will also be distributed between these 
sub-levels. Thus, regardless of strain, the fine structure of the 3 E triplet will average 
to a single splitting of 3D2,a 1 = 1-42 GHz between the m s = and m s = ±1 spin sub- 
levels arising from spin-spin interaction (as depicted in figure [6]). Therefore, as observed 
[T9] . the fine structure of the 3 E triplet at room temperature does not vary between 
NV~ centres and appears as an effective orbital singlet split by spin-spin, similar to the 
ground 3 A 2 triplet. 

The matrix representations of orbital operators in the 1 E and 1 E' singlets fl20l) 
indicate that these states should also be susceptible to the dynamic Jahn- Teller effect, 
although this has not yet been observed. The Jahn- Teller effect is likely to have 
significant consequences for the non-radiative dynamics of these states and influence 
their possible participation in the centre's spin polarization process. Clearly, the model 
of the Jahn- Teller effect in the NV - centre needs to be developed in greater detail to 
fully understand the temperature dependence of the centre's properties and the role of 
the Jahn- Teller effect in spin polarization. 

6. Conclusion 

In this article, the molecular model of the NV - centre's electronic structure was 
developed in full. Through the explicit analysis of the Coulomb repulsion interaction, 
an insight was gained into the many competing factors that influence the energetic 
ordering of the singlet states of the ground and first excited MO configurations. Using 
a proof of the positive definite nature of the Coulomb repulsion separation of the states 
of the ground MO configuration, the simplest consistent orbital structure was adopted 
and the avenues of further investigation and routes to confirmation identified. Having 
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E, 



±1 / 



±1 



±1 EL 



±1 
i_ ±1 







Figure 6. Schematics of the fine structure of the 3 E triplet at low temperature low 
strain (left), low temperature high strain (right), and ambient temperature (centre). 
The solid arrows indicate allowed phonon transitions and the dashed arrow indicates 
the transition allowed once E>i,e.\ spin-spin coupling within the two orbital branches 
in the high strain case is introduced. The plot demonstrates that if there were equal 
population in each of the phonon coupled states, the fine structure would average to 
the central structure with a single spin-spin splitting of 3D 2 ,A! = 1.42 GHz regardless 
of strain. 



obtained the orbital structure, spin-orbit and spin-spin interactions were treated using 
perturbation theory and the centre's fine structure and first order state couplings were 
determined. The results were expressed in terms of the simplest set of one- and two- 
electron reduced density matrix elements, allowing direct insight into which parameters 
were known and how to pursue the remaining unknown parameters. The calculation 
of the matrix representations of the centre's interactions with electric, magnetic and 
strain fields allowed the efficient discussion of the effects of strain on the centre's fine 
structure. The matrix representations also formed the basis of a short discussion of the 
temperature dependence of the centre's fine structure and motivated the central role 
that the dynamic Jahn- Teller effect plays in the centre's dynamics. 
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Table Al. The configuration states expressed in terms of holes rather than electrons. 
Note that the hole formulism has not been used in this article and these expressions 
are provided to assist the readers who have previously used the hole formulism. 
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Table A2. The electronic coupling coefficients correct to first order in spin-orbit and 
spin-spin. Parameters are as denned in table 3J 

(0) (1) (2) 

81,1 = 1 " 

si. 4 = - —2i 



S2,5 



/o^l+f 1,B,2 

52.2 = 1 - 

5 2 .3 = - - iN k k^- 



D E; i 

AT ^-L~ D\ E.2 /OAT m Dl,E,l 

S 2,6 - " N V^E^l -V 2iV ^-E^ 

' a E';0 

S3,3 = N K 

S 3 ,2 = - - lN K K-±^ 



S3,5 = - -iV2N K N, 



X± 



E E 



^ Eea—Ee 



■At 



S3,6 = " 

S3,9 = " ~N K K 

S±A = 1 

S4 1 = - — 2i 



UV K 1\ V Ee-.I-Ee-.O 



Ea i; o 



s 4 ,8 = - iy/2 

S5,5 = N, 



/oat DlE,l AT „A_L— DlE,2 



_ Ee.I 

*5,3 = - -iV2N K N v E Jl EE . - 

•56,6 = N v 

AT ^-L ~~ D\ e,2 /7T \T 



S6,3 = ~ 

S6,5 = - N n rj 
S6Q = - -iN^Nr, 



A AT AT «A||W2r)A_L 
-UV K W V Ee;1-E B . 



57.7 = 1 

58.8 = 1 
S8,l = " 



V2 



"^"VE B ,. -E E ..: 



^±+E>l,E,2 

Eea 



S 8,4 = " ^ E^-Ea^o 

S9,9 = N K 

s 9 ,2 = - -iN K ^- 



a E';0 

sg,3 = - N k k 

S9,5 = - 

s 9.e = - -iN K N v EE ,*l EEil 



11V K 1\ V E E ,. -E E ;i 



The electronic solution of the NV centre 



22 



Table A3. Combined matrix representation of Al and Oa 2 orbital operators 
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Table A4. Combined matrix representation of Oe, x and Oe, v orbital operators 
in the basis . *% E<X , 'Kj ,,- 'Kj ., - 'Kr „- *%1 ■ 'Kj:.,- 'K,: ,- 'K,: ,- 'Kj:..,- 

®7°a 2 ' *9°b,xi ®9°E,y}- To obtain the matrix representation of Oe, x , set O a , x -> 

^(ai||V^||e), O a , y -> 0, Of,^ -> -^(e| |Ve | |e), and 6 ,j/ -> 0. To obtain the matrix 

representation of Oe, v , set O a , y -> ^75 (ai | |Ve| je), O a , x -> 0, 6iJ/ -> ^(e|| VeUe), and 
Ob.x -> 0. 
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Table A5. Combined matrix representation of the com- 

ponents of the total spin operator S in the basis 
/<B so <B so cE> so <B so cE> so <B so <B so <b so (b so <B so (b so (b so 

^9°e k' ^9°e 2/}- ^o obtain the matrix representation of S^, set S x h and 5 a , 
-4- 0. To obtain the matrix representation of S y , set S y — > ft and S" x , S* z — »■ 0. To 
obtain the matrix representation of S z , set S z — > ft, and 5 X , 5 y — > 0. 
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